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The global medial structure of regions in 

James Damon 

For compact regions il in with generic smooth boundary B, we consider 
geometric properties of il wiiich lie midway between their topology and geometry 
and can be summarized by the term "geometric complexity". The "geometric 
complexity" of f2 is captured by its Blum medial axis M, which is a Whitney 
stratified set whose local structure at each point is given by specific standard local 
types. 

We classify the geometric complexity by giving a structure theorem for the Blum 
medial axis M. We do so by first giving an algorithm for decomposing M using 
the local types into "irreducible components" and then representing each medial 
component as obtained by attaching surfaces with boundaries to 4-valent graphs. 
The two stages are described by a two level extended graph structure. The top level 
describes a simplified form of the attaching of the irreducible medial components 
to each other, and the second level extended graph structure for each irreducible 
component specifies how to construct the component. 

We further use the data associated to the extended graph structures to express 
topological invariants of fi such as the homology and fundamental group in terms 
of the singular invariants of M defined using the local standard types and the 
extended graph structures. Using the classification, we characterize contractible 
regions in terms of the extended graph structures and the associated data. 

57N80; 68U05, 53A05, 55P55 



Introduction 

We consider a compact region $7 C with generic smooth boundary B . We are 
interested in geometric properties of i7 and B which lie midway between their topology 
and geometry and can be summarized by the term "geometric complexity". Our goal in 
this paper is first to give a structure theorem for the "geometric complexity" of such 
regions. Second, we directly relate this structure to the topology of the region and 
deduce how the topology places restrictions on the geometric complexity and how the 
structure capturing the geometric complexity determines the topology of the region. 

To explain what we mean by geometric complexity, we first consider . If B is 
a simple closed curve as in Figure 1 (a), then by the Jordan Curve and Schoenflies 
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Theorems, $7 in Figure 1 (b) is topologically a 2-disk, and hence contractible. Hence, 
the region in Figure 1 (b) is topologically simple; however, it has considerable "geometric 
complexity". It is this geometric complexity which is important for understanding shape 
features for both 2 and 3 dimensional objects in a number of areas such as computer 
and medical imaging, biology, etc. 



Figure 1 : Simple closed curve in (a) bounding a contractible region in M? in (b) 

This complexity is not captured by traditional geometrical invariants such as the local 
curvature of B nor by global geometric invariants given by integrals over B or Q. 
Nor do traditional results such as the Riemann mapping theorem, which provides the 
existence of a conformal diffeomorphism of with Vt , provide a comparison of the 
geometric complexity of with . For 2-dimensional contractible regions, a theorem 
of Grayson [14], building on the combined work of Gage and Hamilton [10, 1 1], shows 
that under curvature flow the boundary B evolves so the fingers of the region shrink 
and the region ultimately simplifies and shrinks to a convex region which contracts 
to a "round point". Throughout this evolution, the boundary of the region remains 
smooth. The order of shrinking and disappearing of subregions provides a model of the 
geometric complexity of the region. Unfortunately, this approach already fails in , 
where the corresponding mean curvature flow may develop singularities, as in the case 
of the "dumbell" surface found by Grayson (also see eg Sethian [18]). 

The structure theorem for the "geometric complexity" which we give is based on the 
global structure of the "Blum medial axis" M. The medial axis is a singular space 
which encodes both the topology and geometry of the region (see Figure 2). It is defined 
in all dimensions and has multiple descriptions including locus of centers of spheres 
in VL which are tangent to B at two or more points (or have a degenerate tangency) 
as in Blum and Nagel [2], the shock set for the eikonal/ "grassfire" flow as in Kimia, 
Tannenbaum and Zucker [15], which is also a geometric flow on i7 , and the Maxwell set 
for the family of distance to the boundary functions as in Mather [17]. It has alternately 
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been called the central set by Yomdin [20], and has as an analogue the cut- locus for 
regions without conjugate points in Riemannian manifolds. 

The multiple descriptions allow its local structure to be explicitly determined for regions 
in M"+^ with generic smooth boundaries: M is an n-dimensional Whitney stratified 
set (by Mather [17]) which is a strong deformation retract of Q. The local structure 
of M is given by a specific list of local models, by Blum and Nagel for n = 1 , by 
Yomdin [20] for n < 3 (where it is called the central set) and Mather [17] for n < 6, 
and is given a precise singularity theoretic geometrical description by Giblin [12] for 
n = 2. Results of Buchner [5] show the cut locus has analogous properties for regions 
without conjugate points (and our structure theorem extends to this more general case). 
Furthermore, by results in [7, 8, 9, 6], we can derive the local, relative, and global 
geometry of both i7 and B in terms of geometric properties defined on M. 




Figure 2: Blum medial axis for a region in M.^ 



For example, for generic compact regions in , M is a 1 -dimensional singular space 
whose singular points either have F-shaped branching or are end points. This defines a 
natural graph structure with vertices for the branch and end points and edges representing 
the curve segments joining these points (see Figure 3). This graph structure encodes 
the geometric complexity of the region. Furthermore, Q is contractible if and only if 
the graph is a tree. Then the tree structure can be used to contract the region (also for 
computer imaging a tree structure is a desirable feature, as trees can be searched in 
polynomial time). A computer scientist Mads Nielsen has asked whether in the case of 
contractible C , the Blum medial axis still has a "tree structure". We answer this 
question as part of the structure theorem for the medial axis M for generic regions. 

While the global structure of the medial axis as a Whitney stratified set with specific 
local singular structure does capture the geometric complexity of 11 , it is insufficient in 
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Figure 3: Blum medial axis for a region in M? with associated graph structure 

this raw form to characterize the geometric complexity or to see directly its relation 
with the topology. This is what the structure theorem accomplishes. 

We give a brief overview of the form of the structure theorem (Theorem 4.1). 
The structure theorem for general regions 

At the top level we decompose M into "irreducible medial components" M, which 
are joined to each other along "fin curves". Geometrically this corresponds to taking 
connected sums of 3 -manifolds with boundaries. A simplified form of the attaching 
is described by a top level directed graph T(M) . We assign vertices representing the 
Mi, and directed edges from M, to Mj for each edge component of M,- attached to My 
along a fin curve. This is an "extended graph" in that there may be more than one edge 
between a pair of vertices or edges from a vertex to itself (see Figure 4). The resulting 
space obtained by attaching the M; as indicated by the edges of T{M) is the simplified 
form M of the original M. It is homotopy equivalent to M and has the same irreducible 
components. From the topological perspective, the decomposition is optimal in that the 
homology and fundamental group of M are decomposed into direct sums, resp. free 
products, of those for the M,-. Also, M can be recovered from M by certain "sliding 
operations along fin curves" according to additional data defined from M. 




Figure 4: Graph structure given by decomposition into irreducible medial components with 
edges indicating the attaching along "fin curves" 
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At the second level, we describe the structure of each irreducible medial component 
Mi by an extended graph A(M;), denoted more simply by A/, as in Figure 5. Here, 
there is an analogy with the resolution graph of an isolated surface singularity, whose 
vertices correspond to complex curves, ie compact orientable real surfaces, (with data 
the genera and self-intersection numbers) and edges indicating transverse intersection. 
Alternately, we can view the vertices as denoting multiply-punctured real surfaces 
which are attached to the set of intersection points as indicated by the edges. For the 
irreducible medial components there is an analogous structure but with several added 
complications. 

The extended graph A, corresponding to an irreducible medial component M, has 
two types of vertices: "5-vertices" which are attached by edges to "F-nodes". The 
5-vertices correspond to the medial sheets Sy, which are (closures of) connected 
components of the set of smooth points of the component M, (which are not in general 
the connected components of the set of smooth points of M). They are compact surfaces 
with boundaries (possibly nonorientable). The F-nodes correspond to connected 
components yij of the "F-network" J^, of M, . The F-network of M, is the collection 
of "F-branch curves" together with vertices which are the "6-junction points" where 6 
sheets of the Blum medial axis come together in a point along 4 F-branch curves (and 
is described by a 4-valent extended graph). Each edge from a vertex for Sy to a node 
for represents the attaching of the sheet Sg along one of its boundary components 
to the component . 




yik yi,n yi, 



Figure 5: Extended graph structure for irreducible medial component M, : ^-vertices • 
representing connected medial sheets Sij joined to F-nodes ■ representing components 3^,/; of 
the "F-network" 

There is associated data attached to the vertices and edges of A; ; namely, for 5-vertices 
the genus, orientability, and number of medial edge curves on Stj, for F-nodes the 
4-valent extended graph associated to 3^;<., and for each edge the attaching data of the 
boundary component of Sij to ytk . 

From this structure we obtain two alternate ways to view M, : either as obtained by 
attaching compact surfaces with boundaries along (some of the) boundary components 
or by an associated CW-decomposition. These allow us to relate the graph structures 
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and their associated data for all of the M, to the topological structure of O . 

Topological invariants of M, and hence VL, are computed by Theorem 3.2 and Theo- 
rem 3.4 in terms of M, and T{M) . In turn, those of each M,- are expressed in terms of the 
extended graph A, , and the associated data. This data includes the singular invariants 
of each M, such as the numbers of 6-junction points, medial edge curves, components 
of 3^, , the genera of 5,y, etc. From this we define a single "algebraic attaching map" 
for Mi. From this data, we compute in Theorem 3.2 the topological invariants of 
such as homology groups and Euler characteristic. We then deduce in Theorem 3.4 
the corresponding invariants for M and hence, VL. Also, in Theorem 7.4 we give a 
presentation of 7ri(M,) in terms of the associated data for A,-. These provide topological 
bounds on the geometric complexity of Q . 

The answer to Nielsen's question turns out in general to be no provided we wish to 
classify the medial axis up to homeomorphism. The tree has to be replaced by an 
Algorithm we give in Section 1 for decomposing the medial axis into the irreducible 
components M, . The data of the algorithm requires more than just a graph showing 
which medial components are attached to which. However, the simplified version M, 
which sacrifices some of the detail of the attachings, is described by an extended graph 
r(M), and for contractible regions T{M) is a tree. Then the special form which the 
structure theorem takes for contractible regions in M? (Theorem 5.2) gives a complete 
characterization of contractible regions by the following conditions: the extended graphs 
r(M) and the A,- are trees, the medial sheets have genus (and so are topological 
2-disks with a finite number of holes) with at most one boundary curve representing a 
medial edge curve, a numerical "Euler relation" is satisfied which involves the basic 
medial invariants, and a fundamental group condition holds (Condition 5.1) which 
involves a space formed from 3^; and A, . 

The author is especially grateful to Mads Nielsen for initially raising the question about 
the tree structure for the contractible case, which led this investigator to these questions 
and results. 
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1 Decomposition into irreducible medial components 
Generic local structure of Blum medial axis 

We consider a region ncM? with generic smooth boundary B and Blum medial axis 
M. Then by Mather [17], M can be viewed as the Maxwell set for the family of distance 
functions on B; hence, for generic ^, it is a 2-dimensional Whitney stratified set. Also, 
the generic local structure has one of the following local forms in Figure 6 (see eg [12], 
where Giblin gives a very explicit geometric description). 




(a) edge (b) Y-branching (c) fin creation point (d) "6-junction" 

Figure 6: Local generic structure for Blum medial axes in M.^ and the associated radial vector 
fields to points of tangency on the boundary 

Then M consists of the following: (i) smooth connected (2-dimensional) strata; 1- 
dimensional strata consisting of (ii) F-junction curves along which three strata meet in a 
F-branching pattern and (iii) edge curves consisting of edge points of M; 0-dimensional 
strata consisting of (iv) fin points and (v) 6-junction points, where six medial sheets 
meet along with 4 K-junction curves. Connected components of K-junction curves end 
either at fin points or 6-junction points; while edge curves only end at fin points. We 
refer to the union of F-junction curves, fin points, and 6-junction points as the initial 
Y— network y. 

In the simplest form we can view the medial axis as formed by attaching the connected 2- 
dimensional smooth strata to the 1 -complex formed from the K-branch and medial edge 
curves and the fin and 6-junction points. In fact, this approach misses a considerable 
amount of global structure for subspaces of M. To identify this larger structure we will 
decompose M into "irreducible medial components" M, by "cutting M along fin curves". 
Then we further decompose the irreducible medial components M,- by representing 
them as obtained by attaching smooth medial sheets to the resulting F-network formed 
from the union of F-junction curves and 6-junction points in M, . From this we will 
ultimately give a CW-decomposition to compute the topological invariants of M. 
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Cutting the medial axis along fin curves 

In order to proceed, we first explain how we cut along fin curves. 

Suppose we have on the Blum medial axis a fin point xq. Then near xq we can 
distinguish the "fin sheet" which is the sheet that contains a medial axis edge curve 
ending at xq . We can begin following the F-branch curve from xq , while keeping track 
of the fin sheet which remains a connected sheet as we move along the curve. Eventually 
one of two things must happen: either we reach a 6-junction point or another fin point. 

First, if the F-branch curve meets a 6-junction point, then because the fin sheet is 
locally connected near the 6-junction point, we can follow the edge of the sheet as it 
continues through the 6-junction point. After the 6-junction point we have identified 
both the corresponding continuation of the F-branch curve, and the sheet. We can do 
this for each 6-junction point it encounters. As M is compact, eventually the F-branch 
curve must meet another fin point. We shall refer to the edge of the sheet from one fin 
point to the other as nfin curve. 

At the end of the fin curve, what was identified as the fin sheet (close to the F-branch 
curve) from the beginning may or may not be the fin sheet for the end point. If the sheet 
is a fin sheet at both ends, then we refer to the fin curve as being "essential", while if 
it is only a fin sheet at one end, then we refer to the fin curve as "inessential", (later 
discussion will explain the reason for these labels). Examples of these are shown in 
Figure 7. 




Figure 7: Two possibilities for fin curves on a medial sheet and the results of cutting along the 
fin sheets: (a) essential fin curve (b) inessential fin curve 

An example of a region containing an inessential fin curve is given in Figure 8 and 
might be called a "Mobius board", a surfboard but with a "Mobius band" twist. 




Figure 8: "Mobius boai^d" with a "inessential fin curve" 
Qeometry & Topology 10 (2006) 
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Then to cut along a fin curve beginning from a fin point xo , we locally disconnect the 
fin sheet from the direction of xq by locally adding a point of closure to the fin sheet for 
each point of closure on the fin curve (as followed from xq). We must take care as a 
point on the fin curve can locally be a point of closure for more than one part of the 
sheet. After this step, locally at an added closure point, the fin sheet is now a surface 
with piecewise smooth boundary; see Figure 7. If the fin curve is essential, then the fin 
sheet is now locally disconnected from the remaining two sheets still attached along the 
F-junction curve. While if the fin curve is inessential, the fin sheet is still attached at 
the other fin point to the remaining two sheets; again see Figure 7 (b). 

Then we can take the two remaining sheets still attached along that fin curve and smooth 
them to form a smooth sheet along the curve, with former 6-junction points on the fin 
curve becoming 7-branch points (for another F-branch curve). 

After having cut along the fin curves as described, the former points on the fin curve 
have become altered as follows: 

(1) Fin points become points of an edge of the fin sheet. 

(2) A former F-branch point becomes a closure point on an edge of a fin sheet. 

(3) A former 6-junction point become a (topological) F-network point. 

(4) A F-branch point on a base sheet becomes a (topological) 2-manifold point of 
that sheet. 

Because of (4) above, the end result depends upon a further distinction for essential fin 
curves. A type-l essential fin curve will be one which only intersects other essential 
fin curves at 6-junction points; otherwise, it shares a segment of F-branch curve with 
another essential fin curve, and it will be type-2 essential fin curve (see eg Figure 9). If 
we cut along a type- 1 essential fin curve, then the fin sheet becomes disconnected from 
the other sheets (at least along the curve) and this does not alter any other essential fin 
curve. If we cut along a type-2 essential fin curve, then it will alter the structure of the 
other essential fin curves sharing a segment of F-branch curve with it. Hence, we can 
prescribe the following algorithm for decomposing the medial axis. 

Algorithm for decomposing medial axis into irreducible components 

(1) Identify all type-l essential fin curves and systematically cut along type-l 
essential fin curves (it does not matter which order we choose as cutting along 
one does not alter the fin properties of another). 

(2) After cutting along all type-l essential fin curves, we may change certain 
inessential fin curves to type-l essential ones. If so return to step (1). 
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(3) There only remain type-2 essential fin curves and inessential fin curves. Choose 
an essential fin curve and cut along it. If a type- 1 essential fin curve is created, 
return to step ( 1 ). Otherwise, repeat this step until no essential fin curves remain. 

(4) When there are no other essential fin curves, choose an inessential fin curve 
which crosses a 6-junction point, and cut it from one side until we cut across 
one 6-junction point. 

(5) Check whether we have created an essential fin curve. If so then we cut along it, 
and repeat the earlier steps (l)-(3). 

(6) If no essential fin curve is created, then we repeat step (3) until there are only 
inessential fin curves which do not cross 6-junction points. 

(7) Finally we can cut each such remaining inessential fin curve, producing part of a 
smooth sheet as in Figure 7 (b). 

(8) The remaining connected pieces are the "irreducible medial components" M, of 
M. 

Remark 1.1 The distinct connected pieces created following steps (1) and (2) are 
intrinsic to M; while those created using steps (3) and (4) are not because choices are 
involved. Which choices are made typically depends on the given situation and the 
importance we subjectively assign to how sheets are attached. 

Example 1.2 In Figure 9 (a), we have a contractible medial axis with a pair of type-2 
essential fin curves. Depending on which essential fin curve we choose, 1-2 or 3-4, 
we choose to cut along in step (3), we obtain either (b) or (c), which leads to different 
attachings (and hence top level graph) for the irreducible medial components. An 
alternate possibility would be to cut each fin sheet along the fin curves and view them 
as being attached partially along the edge of a fourth sheet. Again, the exact geometric 
form of M may suggest one choice being preferred over the others. 

Example 1.3 In Figure 10 (a), we have a contractible medial axis with 10 fin points 
1-10, and all fin curves are inessential. Depending on how we choose cuts in step (4) of 
the algorithm, we can end up with 1 , 2, or 3 irreducible medial components. 

If we cut from 4 through the first 6-junction point, then 3-6 becomes an essential fin 
curve, and we cut away the fin sheet M\ as in Figure 10 (b). Then further cutting from 
8 through the first 6-junction point, we create another essential fin curve 2-9. Cutting 
along it creates a second fin sheet M2. The remaining inessential fin curves 1-4, 5-7, 
and 8-10 can be contracted to points on edges of the third sheet M3 . Each of these 3 
medial sheets are then irreducible components. 
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Figure 9: Nonuniqueness of medial decomposition resulting from type-2 essential fin curves: 
(a) is a contractible medial axis with only type-2 essential fin curves, and (b) and (c) illustrate 
the results from cutting along the fin curves 1-2 or 3-4. 



Alternatively, after the first cut, we could have instead cut from 9 , and then from 4 
again, and then only inessential fin curves remain without 6-junction points, so they 
contract to a second sheet, and we only obtain two irreducible components. Thirdly, 
we could have begun cutting from 7, then 8, and then 4 twice and we would obtain 
only a single medial sheet with inessential fin curves, leading to a single irreducible 
component. 




Figure 10: Nonuniqueness of medial decomposition resulting from inessential fin curves: (a) 
is a contractible medial axis with only inessential fin curves, and (b) illustrates the cutting of 
irreducible medial component Mi after first cutting from fin point 4 . 
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Constructing the medial axis by attaching irreducible medial components 
along fin curves 

To reverse the algorithm and reconstruct M from the M,- requires the following additional 
attaching data: 

(1) the list of which segments of medial edge curves which will be attached as either 
essential or inessential fin curves 

(2) the order in which the sheets will be attached 

(3) the curve at the j-th stage to which the attaching will be made as a fin curve to 
obtain the {j + l)-st stage. 

The curves in step (3) may cross multiple components, and need only be well- defined 
up to isotopy; however, this isotopy is for the space obtained at the 7-th stage, whose 
singular set must be preserved by the isotopy. Hence for example, in Figure 1 1 , even 
though the medial axis is contractible, if the attaching curve were on the same side of 
the two holes then the resulting spaces would not be homeomorphic. 



Figure 11: (a) A contractible medial axis, with attaching of component M4 to multiple 
components Mi , M2 , M3 (b) Irreducible medial components for (a) 



2 Extended graph structure and medial decomposition 

We next turn to a simplified version of M using the same irreducible medial components, 
but with simplified attaching. The resulting simplified version of M, which we denote 
by M, will still be homotopy equivalent to M. To define M, we consider an alternate 
way to understand the topological effects of attaching along a fin curve 7 . Instead, we 
isotope the attaching map along the support of the original 7 , by sliding and shrinking 
so it now becomes a fin curve 7' which no longer passes through any 6-junction points 
(which have disappeared after the isotopy as in Figure 12). 

We can achieve this by replacing M by a subcomplex which is a strong deformation 
retract, and only differs from M in a small neighborhood of the fin curve. 
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Figure 12: Isotoping a fin curve so it misses 6-junction points 



Lemma 2.1 (Isotopy lemma for fin curves) Let be a fin curve in M. Then there is 
a Whitney stratified set M' C M, which is a strong deformation retract of M and which 
only differs from M in a given neighborhood of 7 . Furthermore, in that neighborhood, 
7 has been replaced by a fin curve 7' which does not meet a 6-junction point. 

Proof We obtain M' as the result of a series of deformation retractions M D M^^^ D 
^(2) :) . . . 3 M^''^ = M' . Here M^^ D M^+'^ corresponds to either moving the fin 
point along a F-branch curve so it is in a neighborhood of a 6-junction point where it 
has the normal form as in Figure 13 (a), or when M^^ already has this form then M^^^^ 
is a deformation retraction across a 6-junction point as in Figure 13 (b). 

In the first case we may construct the deformation because M is analytically trivial 
along a F-branch curve and so analytically a product in a neighborhood of a compact 
segment of a F-branch curve. In the second case, we may use the normal form for 
6-junction points to deform along the shown region. □ 




Figure 13: (a) Deformation retraction along a F-branch curve. (b) Deformation retraction 
across a 6-junction point. In both cases the deformation is across the dark region. 

We view this process as sliding and shrinking the fin curve along its support, across 
6-junction points until the fin curve no longer crosses 6-junction points. Also, after 
isotoping the fin curves, the resulting spaces are still homotopy equivalent. 

After first isotoping the fin curves, we can more clearly see the topological effect of 
cutting along fin curves. 
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Lemma 2.2 In the preceding situation, let 7' be a fin curve which does not meet a 
6-junction point. There are two possibilities. 

(1) Suppose the hn curve 7' is essential, with hn sheet Sj at both ends. Then when we 
cut Si along 7' , the sheet locally becomes disconnected from the remaining sheets of 
the fin curve as in Figure 7 (a). 

(2) If instead the fin curve 7' is inessential, then we cut all three sheets along the 
fin curve; we obtain a single sheet with boundary edge containing the fin curve as in 
Figure 7(b). This may be alternately obtained by shrinking the fin curve to a point. 

Proof In the essential case, by the triviality of M along F-branch sheets, the fin sheet 
remains a single sheet along the fin curve. When we cut it we disconnect it from the 
base sheet as in Figure 7 (a). By contrast, in the inessential case when we cut the fin 
sheet, the edge curve continues along the cut fin curve, joining up with another edge 
curve again. Hence, the fin points become edge points for a single edge curve. 

As an alternate way to view the process, the fin sheet for one fin point becomes part of 
the base sheet. Then we may follow a path around the second fin point to end up on 
the base sheet for the first fin point. Then we may continue the curve around the first 
fin point to be on the sheet which becomes the fin sheet for the second fin point. If we 
extend these curves so they intersect the edge curve, then together they form a closed 
curve, so that when the fin curve is contracted to a point we obtain a 2-disk, one edge 
of which is the edge curve. □ 

Remark 2.3 This Lemma contrasts cutting along essential versus inessential fin curves. 
When we cut along the essential fin curve we disconnect the fin sheet from the base 
sheet, causing a change in homotopy type as in Figure 7 (a). By contrast, for inessential 
fin curves, when we contract the fin curve to a point p, the sheet then has an edge curve 
which passes through p. In this case, the "fin curve" can be eliminated, and the pair of 
fin points cancelled, without any change in the homotopy type of M. 

Now beginning with a Blum medial axis, we can repeat the Algorithm given in Section 1, 
except at each step, instead of cutting along the fin curve, by the Isotopy Lemma for Fin 
Curves, we slide the fin sheet along the fin curve. We pass all of the 6-junction points 
(which then become F-branch points), until the fin curve lies on a single smooth sheet. 
Then by Lemma 2.2 (b), we may contract the inessential fin curves eliminating their fin 
points. The resulting space with these simplified attachings is M. Finally we cut the 
essential fin curves using Lemma 2.2 (b). The resulting connected components M,- are 
again the irreducible medial components of M. 
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The initial F-network has been altered by the removal of the fin curves to yield the 
F-network y. Then 3^ = (J;-^'' where each 3^,- is the resulting F-network for M,-. 
Then M is obtained from the irreducible medial components by attaching their edges 
which came from essential fin curves to the isotoped positions in some My . There is 
some ambiguity in this construction in addition to that coming from the algorithm. We 
make a choice of the medial sheet that the fin curve passes through. However, the 
irreducible medial components and the homotopy type of M will remain the same. 

This gives rise to a "top level directed extended graph" T{M) whose vertices correspond 
to the Mi , with an edge going from M, to Mj for each segment of an edge curve of M, 
attached to My along a fin curve. 

There is the following relation between the topology of the full medial axis M and that 
of the irreducible medial components {M,}J"^j and the top level extended graph r(M). 

Proposition 2.4 In the preceding situation 

( 1) for any coefficient group G, 

Hj(M; G) ~ Hj{T{M); G) © (©^'^j //y(M,-; G) for all j > 0; 

(2) TTiCM) ~ 7ri(r(M)) * (*^^i7ri(Mi)). 

Of course the graph r(M) only contributes to homology in dimension 1 . Although it is 
an extended graph, we can compute its fundamental group just as for graphs. To prove 
this we also need the next proposition. 

Proposition 2.5 If T is a nonempty connected extended graph, then it contains a 
maximal tree T . Then vri (F) is a free group with one generator for each edge (including 
loops) of r which do not belong to T. 

Proof of Proposition 2.5 First, we construct a connected graph T' from T . For each 
pair of vertices of T joined by an edge we remove all but one of the edges; as well we 
remove any loops (edges from a vertex to itself). What remains is a connected graph V . 
Then such a graph has a maximal connected tree T (by eg Spanier [19, Chapter 2]). 
This is also a maximal connected tree for T . 

There is only one special case which occurs for F-networks, and hence which we must 
allow. It is a single without any vertex. Then we have to artificially introduce a 
vertex so becomes a loop on the vertex. Then the maximal tree is just the introduced 
vertex. However, we emphasize that we will not count such "artificial vertices" for later 
numerical relations involving the number of vertices. 
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Then by repeating the same proof in [19, Chapter 2] as for the fundamental groups of 
graphs, vri(r) is a free group with one generator for each edge (including loops) of T 
which do not belong to T. □ 

We obtain as an immediate corollary. 

Corollary 2.6 If M is contractible, then so is each Mi contractible, and furthermore, 
T{M) is a tree. 

Proof of Corollary 2.6 If M is contractible, then M is connected, tti (M) = , and 

Hj{M) = for all j > 0. Thus, by (2) of Proposition 2.4, 7ri(r(M)) = 0; thus, r(M) 
is a tree. In addition, for each /, by the combination of (1) and (2) of Proposition 2.4, 
7ri(M,) = 0, and Hj(Mi) = for all j > 0. Also, by definition, each M,- is connected. 
Then by the Hurewicz Theorem, TTj(Mi) = for all j > 0. However, is a 
CW-complex, so by another Theorem of Hurewicz, M, is contractible. □ 

Proof of Proposition 2.4 By our earlier discussion, M is homotopy equivalent to the 
space M obtained by attaching M,- to Mj along an edge segment of M,- to a smooth 
sheet of Mj . However, this edge segment can be homotoped to a point so M is instead 
homotopy equivalent to the space obtained by attaching an external curve segment Q,y 
from Xij on the edge of Mj to a point ytj in a smooth sheet of Mj . Then from one fixed 
point x,o on M,- we can choose disjoint curves from x,o to the other points Xy. Then 
we choose a maximal tree T in r(M) and choose the attaching curves a,y corresponding 
to the edges of J, as well as the curves /3,y in each M,-, which give a tree from each 
Xio . The resulting space X which is the union of these trees and the M,- is homotopy 
equivalent to the pointed union of the M, . Hence, 

7ri(X) ~ and HjiX) ~ eUi^jiMi) for j > 0. 

Then as in the computation of tti of a graph as in [19, Chapter 2], we inductively add 
one of the remaining edges and show using Seifert-Van Kampen that we are taking a 
free product with a free group on one generator. After k steps, where k is the number of 
additional edges (including loops) of r(M) not in T, we obtain that 7r\{M) is isomorphic 
to the free product of tti(X) with a free group on k generators, which is exactly the 
fundamental group of r(M) . 

For the case of homology we instead use Mayer- Vietoris in a similar argument. □ 

Proposition 2.4 is valid for any region Q . Thus, we have a relation between the topology 
of M, and hence Q (as M is a strong deformation retract of Q), and the individual 
irreducible medial components M, and the top level extended graph T{M) . In the special 
case that Q is contractible, then so also is M and each M,-, and r(M) is a directed tree. 
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Example 2.7 (Knot complement regions) The simplest example of contractible O 
corresponds to boundary B ^ . A slightly more complicated example would have B 
consisting of two components and a torus . In the case that belongs to the 
compact region bounded by in , then 17 is the region inside and outside of the 
torus . This is a "knot complement region". 

If C is a knot, then the complement C = M?\K is not a region in our sense. 
However, we can take a small tube around the knot and a large sphere S\ containing 
the tube. Then the region i7 = D\\mt{T^) is a compact region with boundary S\ U dT^. 
Also, $1 C C is a strong deformation retract. For a generic knot with e > and R 
sufficiently large, the Blum medial axis of 17 will be generic. 

To see what form the irreducible medial components {M,} take for such a knot 
complement region, we note a few results from knot theory. First, by Alexander duality, 
Hi{Q) ~ Z for / < 3 and otherwise. Also, tti(C) does not split as a free product 
of nontrivial groups. Hence, by Proposition 2.4 there are just three possibilities: (i) 
7ri(r(M)) 7^ so all M,- are simply connected; or r(M) is a tree and all but one of the 
Mi, say Ml , are simply connected; and then either (ii) H2{M\) ~ Z, or (iii) H2{M{) = 
and for a single other Mi , say M2 , H2{M2) ~ Z. 

In the first case, we must have 7ri(C) ~ 7ri(r(M)) ~ Z, which implies K is the unknot. 
In this case, all but one of M, , say M\, have trivial reduced homology and hence 
are contractible, while M\ is simply connected with reduced homology except for 
HiiMi) ~ Z. Then T{M) is homotopy equivalent to , M\ is homotopy equivalent to 
, and M (and hence C) is homotopy equivalent to the pointed union V (which is 
homotopy equivalent to a standard homotopy model for the complement of the unknot, 
the torus with a disk attached to one of the generators of vri ). 

In the second case, T{M) is a tree and all of the other M, except M\ have trivial 
fundamental group and reduced homology, and hence are contractible. There is then 
exactly one noncontractible irreducible medial component M\ , and C is homotopy 
equivalent to it. 

In the third case, all but M\ and M2 are contractible, M2 is homotopy equivalent to , 
and T{M) is a tree. Then C is homotopy equivalent to M\ V M2 ~ Mi V 5^ . 

In all of these cases, all but one or two M,- are contractible and so contribute nothing to 
the topology of the region but do contribute to its geometric complexity. Next, we shall 
determine the structure of the and relate this to their topology. 
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3 Structure of irreducible medial components 

Next, we decompose each irreducible medial component M, , which is connected and 
locally has only F-junction points, edge points, and 6-junction points (but no fin points). 
We may apply the local closure procedure used in Section 1 to cut the sheets of M, 
along the remaining F-network as in Figure 14. We do this by locally adding closure 
points to each of the three sheets meeting at a point of the F-network, or to each of the 6 
sheets meeting at 6-junction points. This disconnects the connected components of the 
set of smooth points of M, from the F-network, turning them into compact connected 
surfaces with piecewise smooth boundaries. We refer to each such compact surface with 
boundary 5 as a medial sheet of M, . The data attached to each medial sheet consists 
of the triple (g, o, e). Here g is the genus of S (after attaching disks to the boundary 
components), o = 1 or corresponding to whether S is orientable or nonorientable (a 
nonorientable surface with at least one boundary component can be embedded in ), 
and e denotes the number of boundary components of S (ie 5^ 's) which are medial 
edge curves of M (keeping in mind that part of the edge might have originally been an 
essential fin curve). 

Then each M, is formed by attaching certain of the boundary components of medial 
sheets in M,- to the F-network J^,- of M,-, reversing the cutting process above. The 
F-network 3^,- has connected components {3^;i}^Li , and we denote the collection of 
medial sheets of M, by {5,y}j!Li . 



Ml 



M — 



M2 




M3 























Figure 14: Cutting the medial axis first along fin curves and then along the F-network y 

We assign an extended graph A, to each such irreducible medial component M, as 
follows. The graph will have two type of vertices: to each sheet Sij will be associated 
an 5-vertex, and to each connected component we associate a vertex which we 
refer to as a F-node. We will associate an edge from the S'-vertex associated to Sij to 
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the F-node corresponding to , for each boundary component of Sjj (ie an 5 ) which 
is attached to . 

Each medial sheet 5,y has the associated data (g,o,e). Second, each connected 
component 3^/<. of the F-network 3^, can be described by an extended graph n,jt. 
Its vertices correspond to the 6-j unction points of yj^, and there are distinct edges 
between vertices which correspond to the distinct F-junction curves of y joining 
the corresponding 6-junction points. Because a 6-junction point has exactly four 
F-junction curves ending at it, each vertex of the graph will have four edges. Thus, 
each vertex will have a valence of 4, and we will refer to such an extended graph as a 
"4-valent extended graph". 




U O .4 



Figure 15: Simplest examples of irreducible medial components, their F-networks and the 
reduced network graphs 

A basic result we need for 4-valent extended graphs is a consequence of the following 
result for ni-valent extended graphs. 

Proposition 3.1 If IT is a nonempty connected m-valent extended graph with k 
vertices, then (m or k must be even and) 7ri(n) is a free group on ^ ■ (m — 2) + I 
generators. 

Hence, for a nonempty connected 4-valent graph 11 with k vertices, vri(n) is a free 
group on ^ + 1 generators. 

Proof An extended graph 11 is still a 1 -complex. Hence, 7ri(n) is a free group on £ 
generators, where £ = rk//i(n). If V, respectively E, denotes the number of vertices, 
respectively edges, of 11, then by the m-valence of 11, 2E = mV. As V = k, 

I = i-x(n) = i-(v-E) 

k 

= \-{k m) 

2 

which equals the value claimed. □ 
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Associated weighted graphs of extended graphs 

In order to work with an extended graph T , we simplify the description by defining 
an associated weighted graph T . T will have the same vertices as T ; however, for 
each pair of vertices with at least one edge between them, we remove all but one of the 
edges. Likewise, we remove all loops. Then we assign to each vertex v,- the number 21 
where there are i loops at v,- . Also for the edge between vertices v, and vj , we assign 
the integer m which is the number of edges in F between v, and vj . Then the valence 
at V, is 2^ + ^ nii where we sum wi, over the edges of T which end at v,- . If for a 
vertex the number of loops is then we suppress it; likewise if there is only 1 edge 
between two vertices, we usually suppress the 1 . For extended graph n,i we refer to 
this graph I!;,;, as the reduced Y-network graph of . For such a reduced graph, at a 
vertex the sum of the vertex number and edge numbers always equals 4 . There are very 
few possibilities for vertices as shown in Figure 16 




Figure 16: The possible local vertex structures (at the enlarged vertex) for reduced F-network 
graphs. The case of 5' is the exceptional case of nonempty F-network without any vertices. 

Finally, we assign this data associated to to the graph A, as follows. The data of Stj 
will be attached to its 5-vertex, the graph 11;^: (or its reduced graph n,yt) assigned to 
the F-node of , and the topological attaching data will be assigned to the edge from 

Sij to yik. 

Next, just as we did for M, we compute the topological invariants of each irreducible 
medial component M, . This time it will be in terms of the extended graph A, and the 
associated data. 

Associated numerical data for irreducible medial components 

The numerical invariants can be divided into those that are invariants of the graph Kj , 
and those which are given by the associated data attached to the vertices and edges of Aj . 
In general, for an invariant Ij of the medial component Mj , we let / = Ij , summed 
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over the iiTeducible medial components of M, be the associated global invariant of 
M. The only exception is our definition of i/{M) given by (3-2). We summarize these 
invariants in Table 1 . 

Graph-theoretic numerical data of Aj 

The key invariants associated to Aj are: the number sj of 5-vertices, the number cj of 
7-nodes and the first Betti number \j of Aj (as a topological space). 

We can express other invariants such as the (reduced) Euler characteristic in terms 
of these. For example, the number of edges of A, = sj + cj + Ay — 1 . We introduce 
another purely graph-theoretic invariant which appears in computing xi^j) '■ 



where the sum is over the irreducible medial components of M and (3i is first Betti 
number of T(M) (as a topological space and = — x(r(M))). 

Numerical associated data of Aj 

Second, we list the numerical data defined from the associated data. First of all, an 
S-vertex is associated to Sjt , which is a compact surface with boundary. The boundary 
will have connected components which are either edge curves or are attached to the 
7-network yj . We let ejk denote the number which are medial edge curves of Mj . For 
each connected component yjk of the F-network, we let vjk denote the number of 
vertices in Iljic. This is the number of 6- -junction points in the F-network yj. 

Next, Sjk has a genus in the sense that if we attach a 2-disk to each boundary component, 
then we obtain a compact surface (without boundary). By the classification of surfaces, 
if it is orientable, it is homeomorphic to a connected sum of tori and the number is its 
genus; while if the surface is nonorientable, then it is a connected sum of real projective 
planes, and that number is the genus. In either case we denote the genus by gj^, and 
refer to it as the genus of Sjt . We define the weighted genus of Sjk to be 



(3-1) 



uiAj) 




Using this we also introduce the global graph-theoretic invariant of M: 



(3-2) 




(3-3) 




2gjk if Sjk is orientable 
gjk if Sjk is nonorientable 
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Invariants of 
medial components 

Graph invariants 

Si 



Global invariants 
of medial axis 



Ay 



Graph data 
invariants 



A 

v{M) 



e 
G 

q 
Q 



Description of medial 
component invariants 

number 5-vertices of A^- = 
number medial sheets of Mj 
number F-nodes of Ay 
= number connected 
components of yj 

y{Mj) given by (3-1); 
v{M) given by (3-2) 



number vertices of lij = 
number 6-j unction points of Mj 
number medial edge curves of Mj 
total weighted genus of Mj (3^) 
sum of qjt given by (3-5) 
first Betti number of associated 
1 -complex Syj C Mj 



Table 1 : Invariants Ij of Medial Components Mj and their Global Versions 
/ for the medial axis, where / = Ij . 



Then the total weighted genus of Mj 

(3-4) Gj =* gjk summed over all medial sheets Sjk of Mj . 

k 

From the weighted genus, we define qju as follows. 

(3-5) a-, - l~SJ' + 'j>'-^ ifejk>0 

Note this is the first Betti number of a 1 -complex in Sjt which will be used in Sections 
6 to 8 to compute the homology and fundamental group of Mj . Then we define 

qj = qjk summed over all medial sheets Sjk in Mj . 
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Although Sj , the number of medial sheets in Mj , is a graph theoretic invariant, there are 
subsets of S'-vertices whose definition depends on associated data. Of these the most 
important is sqj, the number of sheets without edge curves. This is a sum SQoj + s^nj, 
where sqoJ, resp. s^nj, denotes the number of orientable, resp. nonorientable, sheets 
without medial edge curves. 

Then by (3-3), (3^) and (3-5), there is the following relation between these invariants: 
(3-6) qj = Gj + ej - (sj - soj) 

Topological invariants of irreducible medial components 

We next describe how to determine the Euler characteristic, homology, and fundamental 
group of Mj. Li Section 8, we introduce a "minimal 1 -complex" Syj constructed from 
the singular data of Mj. It consists of three contributions: the graph structure Aj, the 
y-network yj, and a 1 -skeleton from the medial sheets. The ikHiiSyj) is given by 
Qj = Xj + {vj + Cj) + qj , which is the sum of the three contributions. Then in Section 8, 
we define an algebraic attaching homomorphism 

From this homomorphism we can compute the homology of Mj in terms of the singular 
data. Second, we shall also define in Section 7 a set of generators and relations to 
compute the fundamental group of Mj . These provide two ingredients for the topology 
ofM,-. 

Theorem 3.2 Suppose that Mj is an irreducible medial component. Then there are the 
following properties in homology: 

(1) Mj has torsion free homology. 

(2) The reduced Euler characteristic of Mj is given by 

XiMj) = soj - Qj. 

Alternately, it can be written in terms of the extended graph Ay, and associated 
data of Aj : 

^3_7^ X(Mj) = i^{Aj)-{Gj + ej + Vj) 

= - ((^j + Vj + Cj + Gj + Ay). 
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(3) The homology groups H2{Mj), respectively Hi{Mj), are the kernel and cokernel 
of the algebraic attaching homomorphism "ifj; see (8-4). Furthermore, there are 
bounds 

rk{H2{Mj; Z)) < so^j and qj < rk{Hi{Mj- Z)) < Qj - sonj. 
In addition, there are the properties for the fundamental group: 

(4) There is a continuous map ipj : My — > Ay such that the induced map 7ri(My) 
7ri(Ay) is surjective. 

(5) The fundamental group has a presentation by generators and relations 

where the relations rji are given by (7-3). 

Proof For (1), it is sufficient to show H^{M) is torsion free; for H^{Mj) is a direct 
summand by Proposition 2.4. Then since H^(M) ~ H^(Qo) where Qo = Q\B, it 
is sufficient to show H^:(Qo) is torsion free. We may decompose M^\i3 = QqU Q' 
as a disjoint union. Thus, //*($lo) is a direct summand of H^(M.'^\B). However, by 
Alexander duality Hj(]R.^\B) ~ H^^-i(B) for j > 0. Since ;S is a compact orientable 
surface, H*{B) is torsion free, giving the conclusion. 

The derivation of the formula for the Euler characteristic, and the computation of the 
homology in terms of the algebraic attaching homomorphism are given in Section 8. 

The computation of the fundamental group is given in Section 7, using the CW-decom- 
position of My in terms of the singular structure given in Section 6. □ 

These results provide the following restrictions on geometric complexity. 

Corollary 3.3 Consider an irreducible medial component My . 
First, suppose H\{Mj) = 0. Then we have the following: 

(1) The graph Aj is a tree. 

(2) Each medial sheet Sjt has genus , ie it is a 2 -disk with a finite number of holes. 

(3) At most one of the boundary components of Sjk is an edge curve of Mj . 

(4) Suppose, in addition, H2{Mj) = 0. Then ej is the number of medial sheets 
with edge curves, and the following "Euler Relation" holds between the number 
of sheets without edge curves (LHS) and the (RHS), which is rk//i(3^y), a 
topological invariant of the Y-network 

(3-8) Sj - ej = Vj + Cj. 
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(5) If Mj is simply connected (without requiring H2{Mj) = 0), ttien Mj will satisfy 
the fundamental group relation given in Condition 5.1. 

(6) Finally, if Mj is contractible, then all of the preceding hold. 

The derivation of the Corollary from Theorem 3.2 will be given in Section 8 

Finally, we combine the results of Theorem 3.2 with the decomposition in Proposition 2.4 
to compute the homology and reduced Euler characteristic of M from the singular 
invariants measuring geometric complexity. We let 

where sq = J2j ^oj is the total number of medial sheets in all of the My without edge 
curves (and Q = T.jQj)- 

Theorem 3.4 Let M be the Blum medial axis of a connected region Q with smooth 
generic boundary B. Then H^{Q) ~ H^{M) have the following properties. 

(1) M (and $7 ) have torsion free homology. 

(2) H2{M; Z) ~ ker(«') and Hi (M; Z) ~ coker {^) © {T(M); Z) . 

(3) ri(//2(M;Z)) < soo and q + (3i < rkiHi{M;Z)) < Q-son + Pi, where 
sqo, resp. SQn denote the number of orientable, resp. nonorientable, medial sheets 
of M without edge curves and as before Pi = rk(Hi(r{M); Z)). 

(4) The reduced Euler characteristic is given by 

X(M) = u{M) - (G + e + v) = s-{e + v + c + G + p) 
where P = Pi 

(5) In the case Q (and hence M) is contractible, G = 0, P = 0, and we have the 
relation 

s — e = c + V. 

Proof The proof follows from M being a strong deformation retract of Q (see eg 
[7]), using the sum formula given in Proposition 2.4, and substituting the formulas for 
Hi(Mj; Z) and xiMj) given in Theorem 3.2. □ 

Example 3.5 The simplest irreducible medial component M,- has 3^,- = . Then M,- is 
a compact connected surface with boundary so Si = 1 and e, = number of boundai^y 
components. The F-network invariants c, = v, = and the graph A; consists of a 
single S'-vertex so A, = . Then G,- = Qi = qi is the rank of Hi (Mi) , and v(Mi) = 1 . 
If Ci = then M, must be orientable, sqi = 1 , and is the zero map; while if e, > 
then soi = and "i>i is again the zero map. 
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Example 3.6 The next simplest case has 3^, = o . Already there arises a question of 
which nonintersecting embedded surfaces can be the medial sheets. The easiest case is 
when the medial sheets are 2-disks and annuli (examples in Figure 2, Figure 14 and the 
first three examples of Figure 15). However, an example exists of a single medial sheet 
which is annulus with a single medial edge curve to form a "collapsed version of the 
umbilic bracelet" as in Figure 17. 

4 General structure theorem 

We can now state the full structure theorem for the Blum medial axis, and in the next 
section state the form it takes for contractible Q . 

Theorem 4.1 (General structure theorem) Suppose i7 C is a compact connected 
region with generic smooth boundary B and Blum medial axis M . Then there is 
associated to $7 a two level graph structure. 

(1) At the top level r(M) is a directed extended graph. It consists of vertices 
corresponding to the irreducible medial components Mi of M and directed edges 
from Mi to Mj corresponding to each connection attaching an edge of Mi to Mj 
along a fin curve. 

(2) To each Mj is associated a graph Aj with two types of vertices: S-vertices ■ 
corresponding to smooth medial sheets Sji^ of Mj , and nodes • corresponding to 
the connected components of the Y-network of Mj. There is an edge from 
Sjk to yji^i for each boundary of Sjk which is attached to yji^' . 

(3) Each medial sheet Sjk is compact surface with boundary (possibly nonorientable). 
Each Y-network component yjj^t can be described as a 4 -valent extended graph 
Hjj^t (or by the corresponding reduced graph ILji^i ). 

(4) At the third level, the extended graph Aj has data assigned to the S-vertices, 
Y-nodes, and edges. To an S -vertex for Sjk is assigned the data (g, o, e) indicating 
the genus g, o = 1 or denoting orientability or nonorientability and e the 
number of Blum edge curves. To a Y-node for yjk' is assigned the 4-valent 
extended graph Ujic/ . To an edge from Sjk to yjt' , the topological attaching data 
of the boundary circle of Sjk to yj^' . 

(5) Furthermore, the graph theoretic data of each Aj, the number of edge curves, 
6-junction points and the weighted total genus satisfy the equation (3-7). 

(6) The fundamental groups of each irreducible medial component are given in terms 
of generators and relations by (7-2), and the homology groups of each irreducible 
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medial component are torsion free and given as t/ie kernel and cokemel of the 
algebraic attaching map (8-4). 

Proof We have established the decomposition into irreducible medial components, 
which in turn can be represented by the attaching of medial sheets to the F-network. 

The computation of the fundamental group in terms of generators and relations will 
be given in Section 7 based on the CW-decomposition given in Section 6, and the 
computation of both the homology and the reduced Euler characteristic will be given in 
Section 8. □ 

Several examples we have already seen illustrate the general structure theorem. 

Example 4.2 In Figure 1 1, M is formed from 4 medial components, of which M2, 
M3, and M4 are topologically 2-disks (3^ = so a = v,- = 0, and Si = e, = 1). If 
we slide M4 onto Mi , then the graph T(M) is a tree. Mi is formed from 5 medial 
sheets, all of genus 0: two are 2-disks, 2 are annuli, and one is a 2-disk with two 
holes. Except for the annuli, each of the other 3 medial sheets of Mi have single medial 
edge curves. The ^-network is empty except for Mi where it consists of 2 disjoint o . 
Ai is a tree, and the attaching of the 2-disks for Mi kill the two free generators for 
TTiCJ^i) U Ai . Hence, Mi is simply- connected. Also, we see z^(Mi) = 5 — 2 — = 3 
and the other z^(M,) = 1 . Hence, xiM) = (3 + 3- l+0)-(0 + 3 + 3- l+0) = 0. 

Thus, H2{M\) = 0, and Mi and hence all of the M,- and M are contractible. 

Example 4.3 The medial axis shown in Figure 18 is irreducible so M = Mi . There 
are five medial sheets (5 = 5) all of which have genus 0: 3 annuli (two of which have a 
medial edge curve) and 2 disks (soG = 0,e = 2,^ = —1 and the F-network consists 
of 2 disjoint o (c = 2, v = 0). Also, Ay is a tree so A = 0. As above, the attaching of 
the 2- disks for Mi kill the two free generators for 7ri(3^i) U Ai so M = Mi is simply 
connected. We see v{M) = 5- 2 + = 3. Thus, xiM) = 3-(0 + 2 + 0)=l. Hence, 
//2(Mi) = Z, and M = Mi is homotopy equivalent to . 

5 Structure theorem for contractible regions 

For contractible regions Q., the structure theorem for the Blum medial axis takes a 
considerably simplified form. We saw in Corollary 3.3 that the medial sheets and 
extended graphs had to have a special form. As well, we define in Section 7 a set of 
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elements r,y in ■K\iSy'^) corresponding to medial sheets without edge curves; see 
(7-3). The computation of the fundamental group in Theorem 7.4 leads to the following 
condition. 

Fundamental Group Condition 5.1 For an irreducible medial component Mi, the 
set of elements 

{rij : all j for which Sij is without edge curves of M, } 
form a set of generators for tti {Syl) . 

Taken together the conditions on the extended graphs being trees, the numerical "Euler 
Relation" and the preceding Condition 5.1 completely characterize contractible regions. 

Theorem 5.2 (Structure theorem for contractible regions) Suppose i7 C is a 
contractible bounded region with generic smooth boundary B and Blum medial axis M. 
Then there is associated to Q a multilevel directed tree structure which determines the 
simplified structure M associated to M . 

(1) At the top level, T{M) is a directed tree consisting of vertices corresponding to 
the irreducible medial components Mi of M, and there are directed edges from 
Mi to Mj corresponding to the attaching of an edge of Mi to Mj along a tin curve 
in M. 

(2) At the second level, to each Mj is associated a directed tree Ay with two types of 
vertices: S -vertices ■ corresponding to smooth medial sheets Sjk of Mj, and 
Y-nodes • corresponding to the connected components of the Y-network of 
Mi. There is an edge from Sjt to yjj^t if a boundary "circle" of Sjk is attached to 
yjk' ■ Given Sjk and yjk' , there is at most one such boundary "circle". 

(3) Each medial sheet Sjk is topologically a 2-disk with a finite number of holes. At 
most one of the boundary circles of a medial sheet is an edge curve of the Mj . 
Each Y-network component yjk' can be described as a 4-valent extended graph 
Yljk' (or by the corresponding reduced graph tljk' . 

(4) At the third level, the tree Ay has data assigned to the S -vertices, Y-nodes, and 
edges. To an S-vertex for Sjk is assigned (h, e) indicating the number of holes 
and Blum edge curves. To a Y-node for yjk' is assigned the 4-valent extended 
graph Iljk' . To an edge from Sjk to yjk' , the topological attaching data of the 
single boundary circle of Sjk to yjk' . 

(5) Furthermore, for an irreducible medial component Mj, the number of medial 
sheets, edge curves, 6-junction points and the total number of connected 
components of the Y-network 3^,- satisfy the Euler relation (3-8), and each Mj 
satisfies the fundamental group relation Condition 5.1. 
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Conversely, suppose we are given a hounded region VL in with smooth generic 
boundary and Blum medial axis M so that: the top level graph r(M) is a tree; for 
each irreducible medial component Mj, the graph Aj is a tree; the medial sheets are 
topologically 2-disks with a finite number of holes, having at most one boundary circle 
an edge curve of Mj ; the numerical invariants of Mj satisfy (3-8); and each Mj satisfies 
Condition 5.1. Then M and (hence) $7 are contractible. 

Proof If Q is contractible, then so is M which is a strong deformation retract of Q. 
Then by Corollary 2.6, T{M) is a tree and each My is contractible. Then by Corollary 3.3, 
each Aj is a tree, and each medial sheet Sjk is topologically a 2-disk (ie has genus 0), 
with a finite number of holes. At most one of the boundary circles of a medial sheet is 
an edge curve of the M. 

Furthermore, xi^j) = 0, yielding the Euler relation (3-8). As 7ri(My) = 0, by 
Theorem 7.4, Condition 5.1 must hold. 

Conversely, by Theorem 7.4, Condition 5.1 implies that iTi{Mj) = for each j. Hence, 
H\{Mj; Z) = 0. As the graphs Aj are trees. Ay = 0, and all sheets have qj = 0. Thus, 
the Euler condition implies by (3-7) that xiMj) = 0. Hence, as HziMj; Z) is torsion 
free, H2{Mj; Z) = 0. Thus, as T{M) is a tree, by Proposition 2.4, M is simply connected 
and H^:{M;Ij) = 0. As M is a CW-complex, a theorem of Hurewicz implies M is 
contractible. □ 

Example 5.3 (Simplest examples) As the simplest examples illustrating the structure 
theorem, we consider those for which the second level trees are Ily = or o . 

(1) "Simple contractible examples" Each Ily = so each Mj is topologically a 
2-disk. This is the simplest type of region. For example, in Figure 9 and Figure 10, the 
medial axis is decomposed into irreducible medial components of this type. In computer 
imaging, the M-rep structure of Pizer and coauthors [1] is based on the region having 
such a simple medial structure. 

(2) The second type of regions would allow both Ily = or o and the only type of 
medial sheets are 2-disks and annuli. Besides the example in Figure 2, we also see the 
example in Figure 14 and the first three examples of Figure 15 are of this type 

(3) Figure 17 is a collapsed version of the "umbilic torus" investigated independently 
by Helaman Ferguson (who has produced sculptures of it) and Christopher Zeeman who 
called it the "umbilic bracelet". The discriminant for real cubic binary forms is a cone 
on this space. The "umbilic torus/bracelet" is obtained by rotating the hypocycloid of 
three cusps about an external axis in its plane, while rotating by . If we collapse the 
hypocycloid onto the three line skeleton as in Figure 17 (a), and rotate then we obtain 
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Figure 17 (b). It consists of an annulus attached to an inner F-j unction curve and an 
edge curve. It has 11^ = o and a single edge curve, and illustrates how only a single 
medial sheet may meet the y-junction curve. 

(4) The next case would have a reduced network graph of the form "• 4" as for the 
fourth example of Figure 15, which is contractible. 

It appears to be a rather difficult question to determine exactly when examples such as 
(3) or (4) can be included as part of larger medial structures with special properties such 
as being contractible. 




Figure 17: (b) Skeleton of an "umbilic torus/bracelet", obtained by rotating the three intersecting 
lines in (a) around an axis while rotating by ^ 

Example 5.4 (Tree structure for a noncontractible region) A region can have tree 
structures at each of the two levels and consist of medial sheets of genus attached as 
prescribed by Theorem 5.2 yet not be contractible. An example is given by Figure 18. 
M is not contractible, and it is only the Euler relation that fails: s = 5 , e = 2, while 

m = 2 and v = Q,^os — e = 'i^2 = v + m. 



6 Representations of the topological structure of the medial 
sheets 

To determine the topology of an irreducible medial component Mj , we give two alternate 
representations of it as either built up by attaching surfaces with boundaries or by 
attaching cells. These give two alternate useful ways to decompose to compute 
homology and the fundamental group. We introduce two subspaces Syj C Sy'j C Mj , 
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Figure 18: Medial axis M of a noncontractible region with tree structures for T{M) and A 

such that Syj will contain both yj and a subspace homotopy equivalent to the extended 
graph Ay. 

To do this we consider a medial sheet Sjk of Mj . We may use the standard representation 
for compact surfaces as quotients of the 2-disk after appropriately identifying edges 
of the boundary as in Massey [16]. In our case, Sjk is a quotient of a 2-disk with 
a finite number of holes obtained by removing the interiors Bi of embedded 2-disks 
Bi,i = 1, . . . £jk + ejk- Here ij^ is the number of the disk boundaries which will be 
attached to yj , and the other ejk will remain edge curves in Mj . We may assume that for 
example, the edges of the boundary of are identified as given by the classification of 
surfaces; see eg Figure 19 (a). 




Figure 19: (a) Sjk identified as the quotient of a 2-disk with a finite number of holes by 
identifying the edges of the boundary; (b) Edges of holes to be attached to y contained in A, 
and deformation retract r,vt . 

Thus, 5;-^ = (D2\(U;B0)/ ~- To be specific, we suppose that the boundaries of 
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Bi,i = 1, . . . Ijk are the ones which will be attached to yj. Up to isotopy, we may 
assume that Bi,i = 1 , . . . £jk are contained in a 2-disk £)' as shown in Figure 19 (b). 
We let Ajii = D'\{[jf^^ fi,). We may construct curves /5, in Aj^ from a point xjo on the 
boundary of Ajk to boundary points xjj of each 6, , so the curves are disjoint from each 
other except at xjq . 

We let Fjk denote the union of the curves /3, and the union of the boundary edges 
dBi, i = I, . . .£jk. Then Tjk is a strong deformation retract of Ajk- We let 

Syj = (yjUiU.Aj,))/ ^ and Tj = (UTy,)/- 

where ~ denotes the equivalence for attaching the interior boundary edges in the Fyyt to 
yj . Since every edge of yj is the image of one of the attached edges of some Sjt , we 
see that yj C Tj C Syj. Also, Fj is a strong deformation retract of Syj by retracting 
each Ajk to Tj^ leaving yj fixed. 

We next let Sjk = Sjk\int(Ajk). 

Then we can attach each Sjk to Syj along Cjk the outer boundary edge of Ajk . Attaching 
all such Sjk recovers Mj . We summarize these statements with the following lemma. 

Proposition 6.1 Let Mj be an irreducible medial component. Then the following hold: 

(1) There is a subspace Tj C Syj of Mj which contains the Y-network yj. 

(2) Tj is a 1 -complex which is a strong deformation retract of Syj . 

(3) If we collapse the connected components of yj in Tj to points, then the resulting 
quotient space is homeomorphic to the second level extended graph Aj. 

(4) There are compact surfaces with boundaries Sjk obtained from the medial sheets 
Sjk by removing Ajk, a 2-disk with holes whose edges are identified to yj. 
Attaching the Sjk to Syj recovers Mj, up to homeomorphism. 

We shall refer to this as the modified medial sheet representation. 
CW-decomposition of M, 

We further simplify the preceding, by introducing an intermediate space Sy'j . 

We divide the medial sheets into two types: those with edge curves and those without. 
If Sjk has an edge curve, then in the earlier notation it has ejk > edge curves (each 
surrounding a hole) in the 2-disk. We denote these by Ejk. We can choose ejk — 1 
disjoint paths Tki from Xko to points on the distinct edge curves for / < ejk . We can 
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then form Eid = Eu U tu . Likewise, we choose a curve r^o from x^q to the boundary of 
, not intersecting the paths T^i . We can form loops aji centered at x^o by following 
Tki, then going around the i-th edge curve counterclockwise, and then following 
backwards back to Xko (see Figure 20). 




Figure 20: The spaces (a) Sjk , (b) Ljk and (c) Ljk 



We let Ljk = Ajk U % U r^o U (^dV ~) 

k<eji: 

and Ljk is defined as Ljk except we replace Ajk by Fjk. Here "~" denotes the 
identification on the outer boundary of the 2 -disk. In the case that 5^^: has no edge 
curves of Mj , then 

Ljk = Ayi- U U (5dV ~). 

We define 

(6-1) Sy- = 53^y U/t where the union is over all medial sheets in My . 

We describe the relation between Mj and Syj in two steps. First, we relate Sjk and Ljk 
(and Ljk ) and then deduce the consequences for Mj and Syj . 

In the case Sjk has no edge curves, Sjk is homeomorphic to the space obtained by 
attaching a 2-disk to Ljk as follows. 

Lemma 6.2 With the preceding notation: 

(1) If Sjk has no edge curves, then Sjk is homeomorphic to the space obtained by 
attaching a 2-disk to Ljk along the path given in (6-2). 

(2) If Sjk has edge curves, then Ljk is a strong deformation retract of Sjk- 

(3) Ljk is a strong deformation retract of Ljk . 

(4) Furthermore, Ljk is homotopy equivalent (rel Ajk) to Ajk V V, (or equivalently 
Ljk is homotopy equivalent (rel Tjk) to Tjk V V/ ), where \/- S^ is a bouquet of 
qjk 5' 's (with qjk defined by equation (3-3)^. 
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Proof If Sjk has no edge curves then we proceed as follows to recover Sjk from Ljk 
by attaching a 2-disk. We cut the annulus Sju = Sjk\mt{/S.jfi) along t^o. We obtain 
a topological 2-disk. Its boundary is mapped to Ljk by first following tuq , then the 
boundary of the 2-disk in a counter clockwise direction, then backwards along t^q 
and finally around the boundary Cjk of Ay,t but in a clockwise direction. We denote the 
counterclockwise path around dD^ by (JjJ^ , and the counterclockwise path around the 
boundary by 5jk- Then a 2-disk is attached along its boundary to Ljk by the path 

(6-2) ijk = Cjk * 5jk where Qu = tuo * Cj,t * m 

followed by the identification (as is common, a denotes the inverse of a path a). 

Next, suppose instead Sjt has ejk > edge curves. Then we cut Sjk = Sjk\viA{IS.jk) 
along the paths t^q and T^i for / < Cjk ■ The resulting space is topologically an annulus 

def 

with one boundary the edge curve E' = Eje-^. of the ej^-th hole. Then we may retract 
the annulus onto Lj^ by retracting along curves from the boundary edge of E' as shown, 
followed by the identification. This gives Lj^ as a strong deformation retract of Sjk . 

Since this fixes Ajk , we may then further retract Ajk onto Tjk , while fixing the rest of 
Ljk . This gives Ljk as a strong deformation retract of Ljk . 

Lastly, the boundary of after identification is a bouquet of S^ 's, where the number 
is 2gk if Sjk is orientable, and gk if 5^^: is nonorientable. 

Finally, returning to Ljk , if we collapse all of the paths r^o and Tti , / < ejk , to the point 
Xjo, then {Ji^g.^^ Eji U Tj U (dD^ j ~) becomes a bouquet of qkj 5^ 's. This collapsing 
induces a projection map from Ljk to Ajk V V^^i which is the desired homotopy 
equivalence. An analogous argument works for the space Ljk . □ 

Since the deformation retraction for Ljk (resp. Ljk) in (3) of Lemma 6.2 fixes Ajk (resp. 
Vjk), we may apply the Lemma to each medial sheet with an edge curve after it is 
attached to yj . If we let Sy'j denote Syj with all of these medial sheets attached, then 
even though each medial sheet Sjk contracts to a bouquet of 5' 's attached to Ajk, as 
syj is path connected, up to homotopy we may bring these bouquets together and 
conclude the following "CW-decomposition of M/'. 

Proposition 6.3 For an irreducible medial sheet Mj , we have the following: 

(1) Mj is homotopy equivalent to a subspace Mj obtained from Syj by attaching for 
each sheet Sjk without edge curve a 2-disk along the closed path given by (6-2). 

(2) syj is homotopy equivalent (rel Syj) to ( Vfii 5') V Syj, and hence to the 
1 -complex ( Vfii S^) \/ Tj (here qj is defined in Section 3). 
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Proof By Lemma 6.2, we can first strong deformation retract eacii siieet Sjk with an 
edge curve onto Ljk . Then we can obtain each sheet Sjt without an edge curve by 
attaching a 2-disk to Ljk. Because this is done while fixing Syj , we may do this all at 
once obtaining the strong deformation retract Mj . 

Also, by Lemma 6.2, each Ljt has as a strong deformation retract Lj^ (rel Ty/t), which 
is homotopy equivalent (rel Fjk) to ( Vf=i V Tj^. Again doing this for each Ljk 
gives ( Vf4i S^) V Tj as a strong deformation retract of Syj. Thus, Mj is homotopy 
equivalent to the space obtained from ( Vf4i '^^) V Tj by attaching one 2-disk for each 
medial sheet Sjt without edge curves. □ 



7 Fundamental group of irreducible medial components 

To compute the fundamental group of an irreducible medial component Mj , we use the 
CW-complex decomposition we have just given in Section 6. We successively use the 
form of the fundamental group for F-network components given in Proposition 3. 1 and 
that of the graph Ay to compute the fundamental group of Syj . We then compute the 
changes obtained by attaching the Ljk to give that for Syj , and lastly determine the 
effect of attaching the 2-disks corresponding to Sjk without edge curves. 

We first compute the fundamental group of Syj as follows. 
Proposition 7.1 With the preceding notation, 

(7-1) miSyj) ~ VrKr,) ~ ^l(Ay)*«=i^i(3^yO) 

where there are r components for yj. 

Proof By the preceding remarks, it is sufficient to compute -KyiJ^j), viewing Tj as an 
extended graph. We choose a maximal tree Tjk for each component yjk and a maximal 
tree Rj for Ay . Furthermore, as the image of an edge dBjj will cover completely each 
edge of yjk which it meets, we may choose the points Xjk so their images lie in Tjk , eg 
they are vertices . 

We also have to treat the special case where one of the yjk is an o without any vertices. 
As in the proof of Proposition 2.5, we choose a point in o as an artificial vertex, so the o 
can be viewed as a loop on the vertex and then the vertex itself becomes the maximal tree. 
As earlier, we do not count such "vertices" in determining vy for numerical relations. 

Then we let tjk denote the union of the curves /3, in Tjk and form 

= (Uf;iU-ry,)/~. 
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Then each tree Tjk is contractible and if we collapse each Tjk to a point then the quotient 
is isomorphic to the extended graph Ay . 

We may lift the edges of the tree Rj to a collection of edges in Fy , by which we mean 
we can choose curves /?,• in the Tjk which correspond to the edges of Rj . We denote the 
union of those curves by R'j . Then we form 

Tj = (/?ju(Ur;,))/~ 

where now "~" only involves attaching the points Xjk of R'j to the corresponding points 
of the Tji. 

Lemma 7.2 Tj is a maximal tree in the extended graph tj . 

Proof First, '7/ is a tree because each Tji is contractible and if we collapse each to a 
point we obtain Rj which is contractible; hence, Tj is contractible and is a graph. 

To see Tj is maximal in Tj , we note any other edge of Tj either corresponds to an edge 
of one of the yj^ , or to another edge of Ay under the collapsing map. In the first case, if 
the edge r/ is in yjk , then as Ty,- is maximal Tji U is no longer simply connected; thus 
nor will be 7y U . If instead the added edge rj maps to another edge rj' of Ay under 
the collapsing map, then Rj U rj' is no longer simply connected. Then we may lift a 
noncontractible loop 7 of Rj U ry' to a loop in 7y U 77 (joining edges of 7 lift to edges 
which share a common F-node, say yjk, and hence can be connected by a path in Tji), 
destroying the simple connectivity of Tj. Thus, Tj is maximal. □ 

Then by Proposition 2.5, 7ri(53^y) is a free group with a free generator for each edge of 
Fy not in Tj. There are Ay edges of Ay, and vyi + 1 edges for each yjk by Proposition 3.1. 
Thus, there are a total of Ay + vy + cy free generators (recall vy = Vjk and cy is the 
number of connected components yjk of 3^y). Such a free group is the free product of 
the free groups on Vjk + 1 generators which are isomorphic to iriiyjk), and a free group 
on Ay generators which is isomorphic to vri (Ay) . □ 

The 7ri(3^y^:) are naturally identified as subgroups of 7ri(53^y). Also, the collapsing map 
(p: Tj ^ Ay is a homotopy equivalence; hence, 7ri(ry) ~ 7ri(Ay). Then the inclusion 
Fy C Syj , identifies a subgroup isomorphic to tt\ (Ay) . 

The second step is to compute the fundamental group of Sy'j . It is homotopy equivalent 
to (V* 1-^0 V53^y. Hence, again by the Seifert-Van Kampen theorem, 

^x{sy'j) ~ T^i{syj)*F^. 

where Fq. is a free group on qj generators corresponding to the generators of Vfii*^^ 
obtained from the Ejk and the identifications on the boundaries dD^ / ~. Hence, 

Corollary 7.3 vri(53^j) is a free group on Qj{= Ay + vy + cy + qj) generators. 
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Computing the fundamental group of an irreducible medial component 

We now use the proceeding to compute the fundamental group of Mj by determining 
the effect of attaching each medial sheet Sjk ■ 




Figure 21: Representing Sjk by attaching a 2 -disk to Syj 

Then we can compute the fundamental group of Mj by attaching the remaining 2-disks 
to Syj ; see Figure 2 1 . First, we recall the loops at x^o : Cjk for each medial sheet without 
edge curves, and 6jk (obtained by following Cjk in the counterclockwise direction). We 
again let Tj denote a maximal tree in Tj and we choose a base point xq £ Tj , and for 
each Xko, a path at in the tree 7y from xq to x^o- As 7^ is a tree, the path homotopy 
class of aic only depends on x^o . We form 

^jk = OLk* Sjk * ak and (jk = ak * Cjk * ak- 

These define elements rjk = 6jk * (j^^ in 7ri(53^j,xo), which after composing with the 
inclusion map yields elements of 7ri(My,xo). 

Theorem 7.4 7ri(My) has the representation 

(7-2) TTiiMj) ~ FqJN 

where Fq. is a free group on Qj generators, and N is the normal subgroup generated 
by the elements rjk = Qk * (Sjk)^^ , one for each medial sheet without edge curves. 
Equivalently, TTi{Mj) has exactly the relations 

(7-3) rjt = 6jk * Cji^^ = I or equivalently 6jk = Qk- 

Also, for Fq. , qj of the generators correspond to the generators from the identified 
boundaries of the 2-disks, and vj + cy generators come from yj. 
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Also, Xj generators come from the generators for vri (Aj) . In fact, the collapsing map 
ipj : Syj — > Aj extends to a continuous map (pj : Mj —>■ Aj which is a surjective map 

on TTl . 

Proof Mj is obtained up to homotopy equivalence by attaching 2 -disks to Syj along 
the closed paths given by (6-1); see Figure 21. By the Seifert-Van Kampen theorem 
this introduces relations killing the corresponding homotopy classes given by (6-1), but 
with the base point moved to jj. These are exactly the relations given in (7-3). The 
identifications of the various generators from the yjk , the edge curves of the Sjk , and 
the identified boundaries of the 2-disks follow from the Seifert-Van Kampen theorem. 

To show that the collapsing map ipj : Syj — > Ay extends to Mj , we first construct 
a quotient Mj of My by collapsing the components 3^y<. to points. Within My is the 
subspace 5J^y, formed from Syj by collapsing the components yj^ to points. The 
collapsing map ipj : Syj Aj factors through Syj . Hence, it is sufficient to extend 
this collapsing map to My and then compose with the quotient map My Mj . 

Then in Syj , Ay^^ has its inner boundary edges collapsed to points to form Ajt , which 
is topologically a 2-disk. As a 2-disk is an absolute retract, we can extend the identity 
map on Ay^; to a continuous map Sjt —>■ Ayjt , where Sjk is Sjk with each boundary 
edge in Ay/; collapsed to a point. Hence, these continuous maps together give a retract 
My —>■ syj . Composed with the collapsing map gives the continuous extension (pj . 

Finally the collapsing map fy —>■ Aj is a homotopy equivalence, and factors through 
(pj : Mj —>■ Aj. Thus, the induced map on (^y* : 7ri(My) —>■ 7ri(Ay) is surjective. □ 

Hence, in the special case when My is simply connected, we obtain Condition 5.1. 



8 Exact sequence for attaching sheets 

Next, we compute the homology and Euler characteristic of an irreducible medial 
component My . We use the two ways of obtaining M : by the medial sheet representation, 
attaching surfaces with boundaries to 53^y, or the CW-representation obtained by 
attaching cells to Syj . Using these we may compute the exact sequences of pairs 
(Mj,Syj) and {Mj,Syj). Furthermore, we have the inclusion of pairs {Mj,Syj) C 
(My, Sy'j) This leads to the homomorphism of exact sequences of pairs using reduced 
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homology. Since all spaces have dimension at most 2 the sequence ends at H2 . 
> H2{Syj) > H2{Mj) > H2{Mj,Syj) 

> H2{Syj) > H2{Mj) > H2{Mj,Syj) 



(8-1) 



HiiSyj) > HiiMj) > Hi{Mj,Syj) 



Hiisy;) > HiiMj) > Hi{Mj,syj) > o 

First, using Proposition 6. 1 and Proposition 6.3 we can compute the groups HiiSyj) 
and Hi{Syj), which are zero for / > 1 and for / = 1 

and HiiSy-) ~ ~ Z^J. 

Second, we can compute the relative groups by excision using Proposition 6. 1 and 
Proposition 6.3: 

(8-2) HiiMj, Syj) ~ ekHi{Sjk,Cjk) 

where the sum is over all medial sheets, and 

Hi{Mj,Syj) ~ ®kHAD\dD^) 

where the second sum is over medial sheets without edge curves. From this, we conclude 
Hi(Mj, Syj) = if j / 2, and 

H2{Mj,sy'j) ~ 

(recall ^o; denotes the number of medial sheets of Mj which are without edge curves). 
Likewise, from (8-2), we compute Hi{Mj, Syj) using the next lemma. 

Lemma 8.1 Suppose N is a compact surface with boundary, with a distinguished 
boundary component C. Let N' denote the surface obtained from N by attaching a 
2-disk to N along the boundary component C . Then 

Hi{N,C) ~ Hi{N'). 

Proof By excision 

Hi{N',D^) ~ Hi{N,C). 
Then by the exact sequence of the pair {N',D^) using reduced homology, we conclude 
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Hence, by Lemma 8.1: 

(1) If has more than one boundary component, or A'^ is nonorientable, then 
H2{N, C) = 0; while if is orientable with only the single boundary component 
C,then//2(A^, C) = Z. 

(2) For (A^, C) = (Syyt, Cjk), rk(Hi{Sjk, Cjk)) = qjk - Sjk, where Sjk = 1 if eju = 
and Sjk is nonorientable, otherwise £jk = 0. 

Thus, by (8-2), Hi{Mj, Syj) = for / > 2, and 

Hi(Mj,syj) ~ I z(*-^-o„,) e (^2^"^- / = i 

where we recall sqoJ , resp. sonj , denotes the number of orientable, resp. nonorientable, 
medial sheets with no edge curves. 

Hence, the homomorphism of exact sequences (8-1) takes the form 
> > H2{Mj) > Z^'O"-' — ^ 

> > H2{Mi) > Z'oj — ^ 

(8-3) 

X(\j+vj+cj) ^ HiiMj) > Z^*--'"".^) © (Z2)-"'">> > 

i 

^(Xj+Vj+cj+qj) , HiiMj) > > 

Recall = © Z'""-^ is the number of medial sheets without medial edge curves, 
and / denotes inclusion of Z*""-' . 

We can draw a number of conclusions from the diagram (8-3). We let Z?,- denote the 
/-th Betti number of My, and we can take the Euler characteristics of either row and 
obtain the same formula. For example, from the first row of (8-3), 

hi - Sooj + (Ay + Vy + Cy) - bi + (qj - So„j) = 0. 

This yields 

XiMj) = soj - (Xj + Vj + Cj + qj) = soj-Qj. 
Substituting (3-6) into this equation yields (3-7) as asserted. 
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Computing homology by the algebraic attaching homomorphism 

Next, we use the second row. The middle homomorphism is exactly the algebraic 
attaching homomorphism 

(8-4) "^j: Z'"J ^ Z^J. 

It sends the generator of H2(Sjk, Cjt) for orientable Sjk without edge curves to the 
image of its boundary (ie Cjt) under the attaching map to yj. Here Z^J represents 
the first homology of Syj . Then by the exactness of the bottom row of (8-3), 
the kernel and cokernel of the "^j are HjiMj), respectively H\{Mj), yielding (2) of 
Theorem 3.2. Furthermore, by the top row of (8-3), we obtain vk{H2{Mf, Z)) < SQoj- 
As H^:(Mj;Ij) is torsion free, we can tensor the top row with Z/2Z and obtain 
rk(Hi(Mj; Z)) > (qj — sonj) + ^Onj = Ij- Also, by the bottom row, or the formula for 
XiMj), 

Tk(Hi(Mj; Z)) = Qj - soj + ^2 < Qj - soj + sqoJ = Qj - sonj. 
This completes the proof of Theorem 3.2 

The third consequence of diagram (8-3) is for the cases when various homologies of 
Mj vanish. 

Proof of Corollary 3.3 First, suppose H\{Mj) = 0. From the first row, we conclude 
that both qj — sonj = and sonj = 0. Hence, qj = sonj = 0. Then qj = qjk is a 
sum of nonnegative integers, so qj^ = 0. By (3-3), each gjk = and if ejt > 0, then 
Cjk = 1 . In particular, Sjk must be a 2-disk with a finite number of holes, and the edge 
of at most one of the holes is an edge curve of Mj . In particular, all of the medial sheets 
are orientable so 

(8-5) sooj = soj = Sj - ej. 

This establishes (2) and (3) of Corollary 3.3. Furthermore, by Theorem 7.4, the 
collapsing map p: Mj ^ Aj induces a surjection on tti and hence on the abelianized 
TTi . Hence, //i(Ay; Z) = 0, so A^ is a tree and = 0, establishing (1) of Corollary 3.3. 

Second, suppose in addition that H2{Mj; Z) = 0. Then 

Thus, SQoj = Ay- + Vj + Cj. 

Hence, as = 0, this equation and (8-5) together imply 

- = V; + 9- 

The remaining conditions on the fundamental group follow from Theorem 7.4. □ 
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9 Algorithm for contracting contractible medial axes 

In the contractible case, the structure suggests a method for contracting a region. We 
would like a specific algorithm for contracting the medial axis when it is contractible. 
First, we can deform it by sliding along fin curves so we have the simplified structure M. 
Then we want to simplify each irreducible medial component. If the medial component 
has a nonempty y-network , then we choose a medial sheet with a medial edge curve, 
and deform the medial sheet so the medial edge curve touches the F-network. At that 
point a transition occurs with two fin points being created as in Figure 22. We can then 
slide the fin curve until it doesn't meet the remaining F-network. This does not change 
the homotopy type of the medial component. 

Figure 22: Creation of a pair of fin points from a medial edge curve 

We repeat this process as long as the F-network of the component is nonempty. Applying 
this to all components, we reach a point where all F-networks are empty. In the process 
we have created more irreducible components attached to each other along fin curves. 
By contractibility, all of the medial sheets are 2-disks. The attaching is described by a 
tree. We then proceed to simplify the tree. 

We take one of the outermost branches. It corresponds to a 2-disk attached to another 
sheet along a segment of its boundary as a fin curve. We contract down the sheet and 
the fin curve to a point. Then we repeat this process until there is only a single medial 
component left. It is a 2-disk and hence contractible. 

The only point about the preceding argument is the assertion that if a medial component 
is contractible then it has an edge curve. This step has not been completed and we are 
left with a conjecture. 

Conjecture If an "irreducible medial component" is contractible then it has an edge 
curve. 

The "irreducible medial component" is in quotes because after the first step we no 
longer know that it really comes from the medial axis of a region. What we really are 
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asking is whether a singular space which is embedded in and has the local singular 
structure of a medial axis satisfies the conjecture. 

10 Discussion and open questions 

There are a number of basic questions which are unanswered. 

(1) The conditions we gave are for a medial axis which already exists. Given an 
abstract model for a medial axis, satisfying all of the properties, this model may 
not be realized in . In particular, we must understand the conditions allowing 
embeddings of the medial sheets so they do not intersect. 

(2) Likewise, there are related restrictions on the attaching maps to F-network so 
that three sheets are attached at each point. What constraints does this place 
on possible structures? This and the preceding point appeal" to involve subtle 
questions in topology. 

(3) In order to complete the algorithm for contracting regions, we need to prove the 
existence of edge curves on contractible irreducible "medial components". 

(4) The decomposition into irreducible medial components requires operations on 
fin curves. These local operations (with the exception of the cutting of essential 
fin curves) are examples of local deformations of medial axes resulting from 
deformations of regions by results of Giblin and Kimia [13] building on the 
results of Bogaevsky [3, 4]. Can we obtain the operations on fin curves as global 
deformations of the region? 

(5) The structure of the medial axis for specific types of regions such as, for example, 
knot complement regions contains topological information about the region. For 
example, if we look for a topological representative of the region for which the 
number of medial sheets is minimal, then as the reduced Euler characteristic 
is fixed, we obtain from (3-7) an expression for sj — xi^j) as a sum of five 
nonnegative integers ej + vj + cj + Gj + Xj . What possible partitions of sj — xiMj) 
into five nonnegative integers are actually possible, and how does this relate to 
the structure of the original knot? 
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